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The Diep and Johnson (DJ) H2-H2 potential energy surface (PES) obtained from the first principles [P. Dicp, 
K. Johnson, J. Chem. Phys. 113, 3480 (2000); 114, 222 (2000)], has been adjusted through appropriate 
rotation of the three-dimensional coordinate system and applied to low-temperature (T < 300 K) HD+o-/p- 
H2 collisions of astrophysical interest. A non- reactive quantum mechanical close-coupling method is used to 
carry out the computation for the total rotational state-to-state cross sections &j 1 j 2 ->j[j' 2 ( e ) an( i corresponding 
thermal rate coefficients kj 1 j 2 _ > .ji y (T). A rather satisfactory agreement has been obtained between our 
results computed with the modified DJ PES and with the newer H4 PES [A.I. Boothroyd, P.G. Martin, W.J. 
Kcogh, M.J. Peterson, J. Chem. Phys. 116, 666 (2002)], which is also applied in this work. A comparative 
study with previous results is presented and discussed. Significant differences have been obtained for few 
specific rotational transitions in the H2/HD molecules between our results and previous calculations. The 
low temperature data for kj^-^j' f {T) calculated in this work can be used in a future application such as a 
new computation of the HD cooling function of primordial gas, which is important in the astrophysics of the 
early Universe. 



PACS numbers: 34.50.-s, 34.20.-b, 34.20.Gj 

I. INTRODUCTION 

The simplest quantum-mechanical scattering prob- 
lem, that between two hydrogen molecules, has been a 
formidable challenge to chemists and physicists working 
in this area. This is because the elastic and inelastic 
H2+H2 scattering have all the complications of a com- 
plex molecular scattering process like rotational and vi- 
brational excitations and de-excitations. Many theoreti- 
cal and experimental methods have been developed and 
used for the study of these processes^— . Additionally, 
there has been great interest in the study of HD+H2 scat- 
tering. Although, this system bears some similarity with 
H2+H2, the identical system symmetry is broken here 
and hence many aspects of the theoretical formulation 
can be tested in the system under different symmetry 
conditions'^— . 

The PESs of H2+H2 and HD-H2 arc basically the same 
six-dimensional functions. The fact follows from a gen- 
eral theoretical point of view and the Born-Oppenheimer 
model2£. However, from an intuitive physical point of 
view, the two collisions should have rather different scat- 
tering output. This is because the H2 and HD molecules 
have fairly different rotational constants and rotational- 
vibrational spectrum. They have different internal sym- 
metries and dipole moments. Collisional properties of 
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these systems are expected to be highly sensitive to dipole 
moments. Further, the HD-H2 PES can be derived from 
H2+H2 by adjusting the coordinate of the center of mass 
of the HD molecule. Once the symmetry is broken in 
H2+H2 by replacing the H by the D atom in one H2 
molecule one can obtain the HD-H2 PES. The new po- 
tential has all parts of the full HD-H2 interaction includ- 
ing HD's dipole moment. However, in a case when a PES 
is formulated for fixed interatomic coordinates between 
the two hydrogen atoms in the H2-H2 system, it would 
be difficult to extract the PES of HD-H2, because the 
position of the HD center of mass is different from that 
of H2. In this work we apply a rotational method for 
calculation of the HD-H 2 PES from that of H 2 -H 2 . The 
method is based on a rotation of the three-dimensional 
space from the body-fixed H2-H2 coordinate system to 
appropriately adjust the HD molecule center of mass. 

The hydrogen molecule plays an important role in 
many areas of astrophysics. For example, the inter- 
stellar medium (ISM) cooling process is associated with 
the energy loss after inelastic collisions of the parti- 
cles. These processes convert the kinetic energy of ISM's 
particles to their internal energies: for instance, in the 
case of molecules to their internal energy of rotational- 
vibrational degrees of freedom. Therefore, in order to ac- 
curately model the thermal balance and kinetics of ISM 
one needs accurate state-to-state cross sections and ther- 
mal rate constants. Theoretical state-resolved treatment 
requires precise PES of H2-H2, and a reliable dynamical 
method for H2+H2/HD collision^^. However, on the 
other hand, experimental measurements of these cross 
sections is a very difficult technical problem. Unfortu- 
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natcly, up to now no reliable experiments are available 
on these collision processes, which have important as- 
trophysical applications. Moreover, different calculations 
with various H2-H2 PESs showed rather different results 
for important H2+H2 thermal rate coefficients. 

The possible importance of the HD cooling in ISM 
was first noted in 1972 by Dalgarno and McCray-2i. Be- 
cause of the special properties of HD this molecule is 
even more effective cooler than H2. This happens at low 
temperatures T < 100 K-2!, where the cooling function 
of HD becomes ~15 times larger than the H2 cooling 
function^ 4 -. As wc mentioned in the previous paragraph, 
to carry out calculation of the cooling function one needs 
to know precise cross sections and thermal rate coeffi- 
cients of the rotational-vibrational energy transfer colli- 
sions. That is why there is a constant interest to reliable 
quantum-mechanical computation of different rotational- 
vibrational atom-molecular energy transfer collision cross 
sections and corresponding thermal rate coefficients 2 ^—. 

A realistic full-dimensional ab initio PES for the H2-H2 
system was first constructed by Schwenke 3 ^ and that po- 
tential was widely used in a variety of methods and com- 
putation techniques. Flower— used Schwcnkc's H2-H2 
PES2£ in a study of H2-HD collision. Later on, new ex- 
tensive studies of the H2-H2 system by Diep and Johnson 
(D J)i and Boothroyd et al. (BMKP)2I have produced re- 
fined PES for the H2-H2 system. These PESs have been 
used in several different calculations 3 - 2 - - — . However, in 
our previous work s 12 ' 3 we found that in the case of low 
energy H2+H2 collision these PESs provide different re- 
sults for some specific state-resolved cross sections. The 
difference may be up to an order of magnitude. This fact 
was also confirmed in work-i 3 -. The BMKP PES, prob- 
ably, needs future improvement s 12 ' 13 ' 34 , because the DJ 
PES gives better agreement with existing experiments for 
H2+H2 than the BMKP potential. 

Nonetheless, as a first trial Sultanov et al. applied the 
BMKP PES to the low-energy HD+H 2 collisions after 
appropriate modificatio n 37 ' 38 . When the results of these 
calculation s 37 ' 38 were compared with prior studies*^ 
a relatively good agreement was found with Schaefer's 
results^. At the same time substantial differences with 
the newer data obtained by Flower was noted^l. There- 
fore, in the light of these circumstances, it would be useful 
to apply another modern H2-H2 PES 1 - to HD-H2 colli- 
sions. The DJ potential was formulated specifically for 
the symmetric H2-H2 system when distances between hy- 
drogen atoms are fixed at a specific equilibrium value in 
each H2 molecule. In this paper we provide the first cal- 
culation describing collisions of rotationally excited H2 
and HD molecules using the DJ PES. We also provide a 
comparison with our previous calculation s 37 ' 38 with the 
BMKP PES, which is a global six-dimensional surface. 

In Sec. |TT] we briefly outline the quantum-mechanical 
close-coupling approach and our method to convert the 
symmetric DJ PES to be appropriate for calculations 
of the HD+H2 system. We represent results for se- 
lected state-to-statc cross sections and thermal rate co- 



efficients for rotational excitation/de-excitation of HD in 
low-energy collisions with 0-/P-H2 in Sec. IIIII Finally, in 
Sec. IIVI we present a brief summary of our findings and 
conclusion. 



II. COMPUTATIONAL METHOD 

A. Dynamical equations 

The Schrodingcr equation for the (12) + (34) collision in 
the center-of-mass frame, where 1, 2, 3 and 4 are atoms 
and (12) and (34) are diatomic molecules modeled by 
linear rigid rotors, is 3 . ' 39 ' 40 
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where is the momentum operator of the kinetic en- 
ergy of collision, R3 is the collision coordinate, whereas 
Ri and i? 2 are relative vectors between atoms in the 
two diatomic molecules as shown in the Fig. [TJ and 
Lft are the quantum- mechanical rotation operators 

of the rigid rotors. Here, M i2 = [(mi + m 2 )(m 3 + 
mi)\/[{m\ + m.2 + mi + TO4)] is the reduced mass of 
the two diatomic molecules (12) and (34) and ^1(2) = 
[m 1 (3)TO 2 (4)]/[( m i(3) + "12(4))] are reduced masses of the 
two molecules. The vectors Ri(2) are the angles of orien- 
tation of rotors (12) and (34), respectively; V(Ri, R2, R3) 
is the PES for the four-atomic system (1234), and E is 
the total energy in the center-of-mass system. The lin- 
ear rigid-rotor model used in this calculation was already 
applied in some previous studie s 21 ' 39 ' 40 . For the consid- 
ered range of kinetic energies of astrophysical interest the 
model is quite justified. 

The eigenfunctions of the operators ^ in Eq. ([1]) 

are simple spherical harmonics Yj iTni (f). To solve this 
equation the following angular-momentum expansion is 
used: 
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where t/„ T „■ t(Rz) are unknown coordinate functions, 
J is total angular momentum quantum number of the 
(1234) system and M is its projection onto the space 
fixed Z axis, and channel expansion functions are the 
following: 
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FIG. 1. (Color online) Four-atomic system (12) + (34) or 
HD+H2, where H is a hydrogen atom and D is deuterium, 
represented by few-body Jacobi coordinates: Ri , R2, and R3 . 
The vector R3 connects the center of masses of the HD and 
H2 molecules and is directed over the axis OZj 9\ is the angle 
between -Ri and R3, 82 is the angle between R? and _R3, ipi is 
the torsional angle, ji , j'2, L are quantum angular momenta 
over the corresponding Jacobi coordinates R\, R2, and .R3. 



here j x + j 2 = ju, 3vi + L = J, fh\, rh 2 , m 12 and m are 
projections of j\, j 2 , ji 2 and L, respectively. The C's 
are the appropriate Clebsch-Gordan coefficients. The 
quantum-mechanical momenta j\, j 2 and L over corre- 
sponding Jacobi vectors R\, R 2 and R 3 are also shown 
in Fig. 1. 

Upon substitution of Eq. ([2]) into Eq. ([T]), one obtains 
a set of coupled second order differential equations for 
the unknown radial functions U^ M (R 3 )^^- 
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where a = (jij 2 ji 2 L). We apply the hybrid modified log- 
derivative- Airy propagator in the general-purpose scat- 
tering program MOLSCAT— to solve the coupled radial 
Eq. flU). 

The log-derivative matrix is propagated to large inter- 
molecular distances R3, since all experimentally observ- 
able quantum information about the collision is contained 
in the asymptotic behavior of functions U^ M (R3 — > 00). 
The numerical results are matched to the known asymp- 
totic solution to derive the scattering S'-matrix S^ a , : 
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where k aa i = \J2M\ 2 (E + E a — E a >) is the channel wave 
number, E a ^ a i^ are rotational channel energies for a and 
a'. The method was used for each partial wave until 
a converged cross section was obtained. It was verified 
that the results are converged with respect to the number 



of partial waves as well as the matching radius, i?3 ma x, 
for all channels included in our calculation. The cross 
sections for rotational excitation and relaxation can be 
obtained directly from the S'-matrix. In particular the 
cross sections for excitation from jij 2 — > j[j 2 summed 
over the final rh! 1 rti! 2 and averaged over the initial m,\rh 2 
are given by 
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The relative kinetic energy of the two molecules in the 
center of mass frame is: 



e = E - BijiCj'i + 1) - B 2 j 2 (j 2 + 1), 



(7) 



where £?i( 2 ) are the rotation constants of rigid rotors (12) 
and (34), respectively, of total angular momentum jim- 
The relationship between the rotational thermal-rate co- 
efficient kj 1 j^^.j' i j! 2 (T) at temperature T and the cor- 
responding cross section 0jij 2 -tj' 1 j' 2 (, e )i can be obtained 
through the following weighted average^ 
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where ks is Boltzman constant and e s is the minimum 
relative kinetic energy of the two molecules for the levels 
ji and j 2 to become accessible. 



B. The BMKP H 2 -H 2 PES 

The BMKP PES 3 -!, is a global six-dimensional poten- 
tial energy surface for two hydrogen molecules. It was 
especially constructed to represent the whole interaction 
region of the chemical reaction dynamics of the four- 
atomic system and to provide an accurate van der Waals 
well. In the six-dimensional conformation space of the 
four-atomic system the PES forms a complicated three- 
dimensional hyper surface 3 ^. To compute the distances 
between the four atoms Ri,R 2 and R3, the BMKP 
PES uses Cartesian coordinates. Therefore it was nec- 
essary to convert spherical coordinates used in the close- 
coupling method^ to the corresponding Cartesian co- 
ordinates and compute the distances between the four 
atoms followed by calculation of the PE S 37 ' 38 . With- 
out the loss of generality the procedure used a specif- 
ically oriented coordinate system OXY ZML. First we 
introduce the Jacobi coordinates {Ri, R 2 , R3} and the 
radius-vectors of all four atoms in the space-fixed coor- 
dinate system OXYZ: {ri, f 2 , f*3, F4} (not shown in Fig. 
QJ. As in Ref. H3 we apply the following procedure: the 
center of mass of the HD molecule is put at the origin of 
the coordinate system OXYZ, and the R3 is directed to 
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center of mass of the H2 molecule along the OZ axis, as 
shown in Fig. [TJ Then R 3 = {R 3 , 3 = 0, $3 = 0}, with 
03 and $3 the polar and azimuthal angles, i?i = r\ — f 2 , 
R-2 = r 4 — r 3 , f\ — £Ri and f 2 = (1 — £)Ri, where 
£ = ma/ (mi + m 2 ). 

Next, without the loss of generality, we can adopt 
the OXYZ system in such a way, that the HD inter- 
atomic vector Ri lies on the XOZ plane. Then the an- 
gle variables of Ri and R 2 are: -Ri = {©i,<I>i = n} and 
R 2 = {02, $2} respectively. One can see, that the Carte- 
sian coordinates of the atoms of the HD molecule are: 

r\ = {xi = £#1 sin 0i, 2/1 =0,zi =£i?icos©i}, (9) 
Ti = {x 2 = -(1 - sin 0i, 2/2 = 0, 

Z 2 = COS©!}. (10) 

Defining ( = 7714/(7713 + 7714), we have 

r 3 = -R 3 - (1-0^2, (11) 
r 4 = R3 + (R2, (12) 
and the corresponding Cartesian coordinates: 

r 3 = {%3 = -(1 - OR2 sin 2 cos $ 2 , 
J/3 = -(1 - O-R2 sin0 2 sin<E>2, 
% = i?3 - (1 - C) R 2 cos© 2 }, 



r 4 = {x A = QR 2 sin 2 cos $ 2 , 
?/4 = C-^2 sin 2 sin <3> 2 , 

Z± = R 3 + C-R2 COS 2 }. 



(13) 
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C. The Modified Diep and Johnson H 2 -H 2 potential 

Below we briefly present our method to convert the 
symmetric DJ H2-H2 PESi to be suitable for the non- 
symmetric system HD+H2. The method is based on a 
mathematical transformation technique, i.e. a geometri- 
cal rotation of the three-dimensional (3D) space and the 
corresponding space-fixed coordinate system OXYZ^, 
as shown in Fig. [2] It is more convenient here to use a 
slightly different from Fig. [1] orientation of the coordi- 
nate system OXYZ: now the center of mass of the H2 
molecule is set at the origin of the space-fixed OXYZ, 
and R 3 (R 3 ) is directed to the center of mass of H2 (HD). 

The few-body system (1234) can be characterized by 
four radius- vectors: {71,72,73,74}, or alternatively by 
so-called three Jacobi vectors: {^1,^2,-^3}- Usually 
the second option is more convenient for describing the 
quantum-mechanical few-body systems. Next, the initial 
geometry of the system is taken in such a way that the 
Jacobi vector R 3 connects the center of masses of the two 
molecules and is directed over the OZ axis. We can also 
choose OXYZ in such a manner that the Jacobi vector 
R 2 lies in the X-Z plane. Finally, the vector R x can be 
directed anywhere. Then the spherical coordinates of the 




FIG. 2. (Color online) The geometrical configuration of the 
HD+H2 system in the framework of two different Cartesian 
coordinate systems: OXYZ and OX'Y'Z'. Here: x is the 
distance between Oh 2 and Ohd, r) is the rotation angle over 
the OZ axis, i.e. the angle between the axes OZ and OZ' 
and between the axes OX and OX' , r*i, 7*2, 73 and F4 are 
the radius-vectors of the corresponding atoms. The vectors 
R2, R3, and R' 3 , belong to the XOZ plain. 



Jacobi vectors are: R\ = (i?i, 6%, ^12), R% = (R 2} 2 ,0), 
and R 3 = {R3, 0, 0). Thus the 4-body system H 2 -H 2 /HD 
can be fully determined with the use of six variables. 
However, the DJ PES has been prepared for the rigid 
monomer model of ^-H^i, so actually we have only four 
active variables in this consideration: R 3l 6\,6 2 , and <pi2- 
The DJ PESi has been prepared as the following func- 
tion for the distance R 3 between two H2 molecules, two 
polar angles #1(2) and one torsional angle 4>i 2 : 



v(R 3 ,e 1 ,0 2 , ( f )12 ) 
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where the angular functions Gi lt i 2t i arei: 
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(15) 



Gooo( 

G202( 
Go22( 
G224( 



'12 



) = 1, 



iia)= -(3cos 2 0!-l), 
H 2 ) = ^(3cos 2 6 2 -l), 
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(3 cos 2 9 2 — 1 ) — 16 sin 6\ cos 9\ sin 8 2 
x cos 6*2 cos 0i2 + sin 2 8\ sin 2 9 2 cos(2cf>i 2 )]. 



(16) 
(17) 

(18) 



(19) 



The coordinate function i 2 i(R 3 ) has been tabulated 
in worki. 
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First we start with the original DJ PES of Eqs. 
(fT5j) — (TH)| . Then we replace one hydrogen atom "H" with 
a deuterium atom "D" . Thus we break the symmetry by 
shifting the center of mass of one H2 molecule to another 
point, that is from Or 2 to Ohd as shown in Fig. [2] 
we also need to take into acount the difference between 
the number of rotational states in H2+H2 and HD+H2. 
The length of the vector x is x — |i?2|/6. Now we ro- 
tate the OXYZ coordinate system around the OY axis 
in such a way that the new OZ' axis goes through the 
point Ohd- The OY' axis and the old OY axis are par- 
allel, and the angle of this small rotation is 77, see Fig. [5] 
This transformation converts the initial Jacobi vectors in 
OXYZ: Ri = {fli.fli.^ia}, R 2 = {i?2,6» 2 ,0} and R 3 = 
{i?3,0,0} to the corresponding Jacobi vectors with new 
coordinates in the new O'X'Y'Z': R[ = {Ri,0' lt (f/ u }, 
R% = {R' 2 ,8' 2 , 0} and R' 3 = {R' s ,0, 0}. As a result of this 
simple procedure we obtain a new PES, namely: 



Vg? (ifcA^a, 0u) -+ vg^Ci^^i,^,^). (20) 

It is quite obvious, that the rotation does not affect the 
coordinate function i 2 ;(i?3) in Eq. (|T5j) . Note, in 
a previous calculation of low energy collision between 
monodeuterated ammonia NH2D and a helium atom, a 
somewhat similar spatial rotational-translational proce- 
dure has been applied to the original PES of the NH3-HC 
system^. 

Next, any rotation of the 3D OXYZ coordinate system 
can be represented by Euler angles, i.e. {a, (3, 7] 45 i 46 . In 
this work we choose the following Euler angles: {a = 
0, P = 77, 7 = 0}. To calculate the value of the rotational 
angle 77 in Fig. [2] one can consider the internal triangle 
AOhdOOh 2 which is shown in Fig. [2J The angle 77 is 
determined from the following formula: 



t?HD 



cot X] 



(R'z + xsm 6' 2 ) 
a; cos 9' 2 



(21) 



The derivation of (|2"Tj) can be expressed as follows. First, 
the angles of the triangle AOhdOh 2 satisfy the follow- 
ing equation (it — 62 ) + 77 + 6*3 = 71 or $2 = r) + 6' 2 . 
Secondly, based on the law of sines for AOhdOOh 2 - 
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where e = ir — 62- Because sin e = cos 82 we have 
x R R' 
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and finally: 



sin rj sin a 2 cos 02 



cos(rj + 0' 2 ) _ R' 
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from which one can directly obtain the expression (|2"TT) . 
In such a way the rotation of the coordinate system from 



OXYZ to O'X'Y'Z' makes a corresponding transforma- 
tion of the coordinates of the atoms in the 4-body system 
and the distance between two molecules. One has the fol- 
lowing relations among old and new variables^-' 
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— sm(9[) sin(?7) cos(0' 12 ), (25) 
_ v ^, _ v „ -sin(02)sin(77)cos(0i 2 ), (26) 

Slllf 77 } 

cot(^ia) = cot(0i) cos(r?) + cot(^) . .\ \ , (27) 
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) = cos(8[) cos(n) 
1) = cos(6> 2 ) cos(77) 

sin(?7) 



R 3 = 



Rf - 2xR' 3 cos(6> 2 ). 



In the calculation of HD+H2 with the DJ PES one has 
to use new coordinates 0' lt 6' 2 , <f>' 12 , R 3 - However, the po- 
tential (fT!))) has been expressed through the old H2-H2 
variables, hence it is to be transformed to new variables 
using Eqs. (|25 l - (|28 l) . 



III. NUMERICAL RESULTS 

Results for the low energy HD+o- /p-H-2 clastic scatter- 
ing cross-sections and few selected quantum-mechanical 
rotational transitions together with the corresponding 
results from previous studie a 37 ! 38 are presented below. 
The cross sections are also compared with the corre- 
sponding results from*. Additionally, we compare our 
results for the thermal rate coefficients © with previous 
calculations**^*. The results for the low energy elastic 
scattering cross sections cannot be compared with other 
theoretical/experimental data. To the best of our knowl- 
edge such calculations do not exist. 



A. Elastic scattering 

In Fig. [3] we show the present elastic scattering cross 
sections for HD+0-/P-H2 collisions at low and ultra-low 
energies calculated using three different potentials: the 
BMKP potential of Sec. HTBl the modified DJ potential 
of Sec. Ill CI and the original DJ potential appropriate 
for the H2-H2 system. In each cross section there is a 
prominent resonance. In Fig. [3] although the shapes of 
three cross sections are similar to each other the reso- 
nance peak and the position of the resonance differ sig- 
nificantly when we use the original DJ potential without 
the modifications described in the Sec. [TTJ This is be- 
cause the original, symmetrical DJ potential* does not 
have all the asymmetrical features of the HD+H 2 interac- 
tion. Moreover, these features are of crucial importance 
for HD+H2 scattering. In Table Hill we show the elastic 
cross sections a e i (e) for a few selected energies calculated 
with the three potentials. We also include in Table lllll thc 
results for scattering lengths do's. As can be seen from 
Fig. [3] and Table IIIII we have obtained extremely good 
agreement between our calculations with the BMKP and 
with the modified DJ PESs. For the considered range 
of the kinetic energies we obtained full numerical con- 
vergence, for instance, the total angular momentum J 
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TABLE I. The elastic scattering cross sections a e i (10~ 16 cm 2 ) at selected relative kinetic energies e (cm -1 ) and corresponding 
scattering lengths a SC att (10 -8 cm) in the o/p-Hb + HD — ¥ 0-/P-K2 + HD low energy collisions calculated with three different 
potentials: modified DJ PES from Sec. Ill CI and original BMKP^i and Dji PESs. Numbers in parentheses are powers of 10. 
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6.0(-l) 


1641.2 


1577.6 


433.8 


1653.2 


1611.5 


437.0 


1.0 


592.0 


520.1 


419.3 


593.3 


522.6 


421.0 


5.0 


504.4 


492.9 


501.5 


504.4 


492.8 


501.8 


10 


331.1 


310.7 


353.3 


331.1 


310.6 


353.1 


100 


87.6 


84.6 


93.2 


87.6 


84.6 


93.1 








Scattering length: a sca 


tt xHT s (cm) 






0.0 


5.27 


5.50 


5.46 


5.70 


5.50 


5.47 



was used up to the maximum value J max = 12 in this 
calculation. 



B. Non-elastic channels 

The main goal of the present study is not to obtain 
results for every possible transition cross section in the 
HD+H2 collision, but rather to demonstrate how the ap- 
propriately undertaken 3D rotation of the symmetrical 
surface of the H2 - H2 system could be adapted for col- 
lision of the H2 and HD molecules. We carry out calcu- 
lations for a fairly wide region of the collision energies, 
i.e. from 3 K to up to 300 K. This temperature interval 
is relevant for future calculations of the important HD- 
cooling functional. Here we compute few transition cross 
sections in which we noticed substantial differences be- 
tween our results obtained with the two different PESs 
from Refs. [l] and[3l| and also between our results and the 
data from previous studies^i^i. 

In Figs. S] to [5] we show these results for a few selected 
state-to-state HD+o-/j?-H 2 integral cross sections. The 
results have been calculated with two different potentials, 
specifically, with the BMKP PE$2± of Sec. MM and the 
modified DJ PES of Sec. Ill CI When possible we com- 
pare our results with existing previous calculations for 
the state-selected total cross sections and thermal rate 
coefficients. It is necessary to mention that the original 
DJ surface^ cannot be correctly used for the non-elastic 
or transition channel calculations, i.e. for the rotational 
state transitions in an HD+H2 collision. If the original, 
unmodified D J potential is applied one can get extremely 
low numbers (~ 10 -35 ) for the HD+H2 rotational state- 
to-state probabilities and cross sections. That is why it 
was necessary to undertake the geometrical modifications 



to the DJ surface described in Sec. Ill CI We obtained a 
fairly good agreement between our calculations with the 
use of the two PESs. Besides some qualitative differences 
in the state-resolved total cross sections the overall be- 
havior of the cross section Vj 1 j 2 -yj' 1 j' 2 ( v ) was found iden- 
tical, where v is the relative velocity of the two molecules. 
However we obtained substantial differences between our 
u 3132^3' 3' ( v ) cross sections and corresponding data from 
Ref . i. 1 2 

In Fig. 2]we plot the total cross sections for the rota- 
tional transitions (02) - (20) and (13) - (11). The notation 
of the rotational quantum numbers can be understood by 
comparison to the following equations: 

HD(0) +H 2 (2) -> HD(2) +H 2 (0), (29) 
HD(1) + H 2 (3)->HD(1) + H 2 (1), (30) 

where the numbers in parenthesis denote rotational quan- 
tum numbers ji, j 2 etc.. On the upper plot it is seen 
that while the two results of this study are in fairly good 
agreement between each other we obtain substantial dis- 
agreements with the result of Ref. [^. The same is true 
in the lower plot, although the behavior of these cross 
sections has a common character. 

In Fig. [5]we show the cross sections for rotational tran- 
sitions in HD+H 2 for the BMKP and the modified DJ 
PESs for the processes HD(1) + H 2 (3) -> HD(0) + H 2 (l), 
and HD(1) + H 2 (3) -> HD(2) + H 2 (l). The correspond- 
ing results from Ref. H are also shown. We again obtain 
a fairly good agreement between our results computed 
with the BMKP and the modified DJ PESs, however, 
differing significantly from the corresponding Schaefer 
result^. Unfortunately we cannot compare these cross 
sections with the results of the calculation by Flower—, 
in which a different H2-H2 potential from 3 ^ was used. 
This is because Flower's data includes results mostly for 
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FIG. 3. (Color online) Elastic scattering total cross sections 
for HD+0-H2 (upper panel) and HD+p-H 2 (lower panel) at 
different kinetic energies e with the BMKR^I potential of Sec. 
HTBl the modified Dji PES of Sec. |nC] and the original DJ 
PES. 



FIG. 4. (Color online) Total cross sections for transition (02) 
(20) and (13) -»■ (11), i.e. HD(0)+H 2 (2) -»■ HD(2)+H 2 (0) 
(upper panel) and HD(1)+H 2 (3) HD(1)+H 2 (1) (lower 
panel) for different velocities v. Present calculations with 
the BMKP£ and modified DJ PESs are compared with those 
of Ref . El. 



the thermal rate coefficients. However, within the next 
subsection in Table II we compare our few selected ro- 
tational state-resolved thermal rate coefficients with the 
corresponding results from Refs. and[2lT 

From the astrophysical point of view, perhaps, 
one needs only precise rotational and in some rare 
cases vibrational state-to-state thermal rate coefficients 
kj 1 j 2 ^j' i j^(T) in H2+H2, HD+H2 etc. These quantities 
are less sensitive to interaction potentials. However, the 
overall behavior of all possible state-selected cross sec- 
tions should be very important for calculation of the ther- 
mal rates as seen in Fig. 2] It appears that the cross 
sections are much more sensitive to the PESs used in the 
calculations. Hence it is useful and even probably impor- 
tant in some specific cases to compare the cross sections 
from various calculations where different PESs have been 
used. 

Further, Figs. [5] and [7] exhibit our results for the state- 
to-state rotational cross sections in transitions HD(2) + 
H 2 (l) -> HD(0) + H 2 (l) and HD(1) + H 2 (2) -+ HD(2) + 
H 2 (0) and in transitions HD(2) +H 2 (0) ->■ HD(0) + H 2 (0) 
and HD(2) + H 2 (0) -> HD(1) + H 2 (0), respectively. We 
obtained fairly good agreement between our own results. 
Additionally, in these rotational transitions fairly good 
agreement with the corresponding cross sections from 
Ref. has also been obtained. Finally, Fig. [5] shows our 



integral cross section for the transition HD(0) +H 2 (2) — ^ 
HD(2) + H 2 (0). This process is also interesting, because 
the transition occurs in the two molecules simultane- 
ously. We would like to name such processes as double- 
transition processes. 

C. Thermal rate coefficients 

In Table II our rotational thermal rate coefficients 
kj 1 j 2 ^j , 1 j 2 (T) are presented. In this paper we choose only 
three rotational double-transitions in the collision: (02)- 
(20), (12)-(20), and (20)-(02). These results were ob- 
tained from corresponding state-resolved integral cross 
sections Cj 1 j 2 ->j' 1 j 2 ( e ) with the use of the expression ([5]). 
In a previous study^i, Flower compared his rotational 
transition thermal rate coefficients with the correspond- 
ing Schaefer data£ and found substantial differences be- 
tween his results and results from Ref. [B| even at high 
temperatures. The point is that in these processes the 
transition occurs in the two molecules simultaneously. In 
such rotational transitions the probabilities and cross sec- 
tions should be very sensitive to the interaction potential. 
Perhaps, because of this reason the results of Refs. d and 
l2ll differ so dramatically for transitions like (02) - (20). 
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TABLE II 

temperatures T (K) in the p-H 2 (a) + HD(/3) 



Results for three selected state-to-state rotational thermal rate coefficients kj 1 j 



P-H 2 (q') + HD(/?') collisions. 



j[j' 2 {T) cm 3 s 1 at various low 

Calculations with different PESs: the original 

21 



BMKP PES*" and the new modified DJ potential from this work. The corresponding older data from other authors*""* are 
also included in this table. Numbers in parentheses are powers of 10. 



T(K) 



Rotational Thermal Rate Coefficients: k 



(T) cm' s" 



(02)->(20) 



(12)^(20) 



(20)->(02) 



BMKP Mod. DJ Ref. 21 Ref. 5 



BMKP Mod. DJ Ref. 21 Ref. 5 



BMKP Mod. DJ Ref. 21 Ref. 5 



10 

20 

30 

50 

70 

100 

200 

300 



1.15(-11) 
9.57(-12) 
9.05(-12) 
8.81(-12) 
8.83(-12) 
8.93(-12) 
9.17(-12) 
9.19(-12) 



6.38(-12) 
5.51(-12) 
5.40(-12) 
5.51(-12) 
5.68(-12) 
5.90(-12) 
6.29(-12) 
6.43(-12) 



7.21(-15) 
6.72(-15) 
7.40(-15) 
9.3(-13) 1.09(-14) 
1.70(-14) 
9.9(-13) 3.21(-14) 
1.1(-12) 1.44(-13) 
3.46(-13) 



1.50(-13) 
1.27(-13) 
1.25(-13) 
1.34(-13) 
1.52(-13) 
1.85(-13) 
3.25(-13) 
4.77(-13) 



9., 



32(-14) 
62(-14) 
53(-14) 

5.64(-14) 7.9(-14) 1.35(-13) 
89(-13) 
86(-13) 
8.78(-14) 4.3(-13) 7.31(-13) 
1.06(-13) 1.28(-12) 



6.35(-14) 
5.43(-14) 
5.34(-14) 



6.06(-14) 
6.72(-14) 



1.6(-13) 



9.. 



1.06(-17) 3.43(-18) 2.53(-20) 

9.18(-14) 4.11(-15) 1.26(-17) 

8.80(-14) 4.48(-14) 1.13(-16) 

5.47(-13) 3.14(-13) 7.5(-14) 8.86(-16) 

1.21(-12) 7.38(-13) 2.82(-15) 

2.23(-12) 1.42(-12) 2.8(-13) 9.15(-15) 

4.58(-12) 3.10(-12) 6.0(-13) 7.69(-14) 

5.78(-12) 4.02(-12) 2.28(-13) 
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FIG. 5. (Color online) Same as Fig. [4] for transitions (13) 
(01) and (13) -> (21). 



FIG. 6. (Color online) Same as Fig. [4] for transitions (21) 
(01) and (12) -> (20). 



In fact, we also obtained substantial differences between 
our thermal rates and with corresponding results from 
Refs. [H and [2l|. Particularly the large differences were 
detected at the lower temperature regime. 

In Table II we show our thermal rate coefficients 
khh^jiti ( T ) for the P rocesses HD(0)+H 2 (2) -> HD(2) + 
H 2 (0), and HD(1) + H 2 (2) -> HD(2) + H 2 (0), and those 
for HD(2) + H 2 (0) -> HD(0) + H 2 (2). We present our re- 
sults calculated with the BMKP— and with the modified 
DJ-^ PESs. As before, our results computed with these 



two potentials are close to each other. However, one can 
see that our results and results from Refs. and Hi] differ 
significantly. This happens particularly at low tempera- 
tures, for instance from 10 K to 50 K. Finally in this sec- 
tion, for the process HD(2)+H 2 (0) -> HD(0) + H 2 (2) the 
difference between our k^^-^j' j^(T) and the results from 
Ref. is about 3 orders of magnitude at T=10 K. The 
reason of this substantial deviation is not clear, although 
it might be a result of using different HD+H 2 potentials 
in the current computation and in Refs. and[2ll where 
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FIG. 7. (Color online) Same as Fig. |4]for transitions (20) 
(00) and (20) -» (10). 
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In the case of the inelastic scattering a + b ^ a' + b' 
the detailed balance principle^! relates the direct a + b 
and the reverse a' + b' processes and their cross sections 

a iah-H'Ji ancl °iU-tf.JV 

(2j a + l)(2j b + \)pl^ b af a]b ^ u ,{E) = (2j' a + 1) 

x(2^ + 1)pL4^«( £ )- ( 31 ) 

Here, E is the total energy, j a (fy and j' a ^ are the initial 
and final rotational quantum numbers, Pa< a ')-^b(b') arc 
the initial and final relative momenta between a and b 
species, Cj^j^fj' ( E ) and a i'j' b ^j a jb^ are dircct and 
reverse integral cross sections respectively. The same 
type of the relationship can be obtained for the thermal 
rate coefficients kj 1 j 2 - > .j' i ji (T). Let us rewrite Eq. (|5J) in 
the terms of the total energy E. Considering that the 
relative kinetic energy e between a and b is Eq. (J7J), we 
compute the total energy from the lowest possible level 
between two channels. If it is associated with the second 
channel a 1 + b', i.e. j^j'^ pair, the formula (JSJ becomes: 

x Pk e ~(E-AE)/k B T dE ^ (32) 
Mia 

Here, AS = \B in (jj + 1) + B 2]2 { ]2 + 1)] - [Br $ + 1) + 
+ 1)] i s the energy gap between the direct and 
reverse channels, and e = p 2 ab /M\ 2 is the kinetic energy. 
Now, for the reverse channel the thermal rate coefficient 
is: 

k fX^h ( T ) = Jj^f / "j^-AA (*0 

x P^V e -E/k B T dE ^ (33) 

Mi 2 

Comparing Eqs. (j3"2"j) and (j3"3")l and taking into account 
Eq. ([3~Tj) we obtain the detailed balance formula for the 
thermal rate coefficients: 



FIG. 8. (Color online) Same as Fig. [4] for transition (20) 
(02). 



the authors also used two different PESs and different 
quantum-mechanical methods. However, we would like 
to point out here that our results computed with two 
newer PESs^i are relatively close to each other. 

D. Application of the detailed balance principle 

By using the time reversibility (reciprocity) principle 
one can obtain the detailed balance equation for the di- 
rect and reverse energy transfer processes or reactions^. 



(2ji + l)(2j 2 + l)kf lj2 ^ fij ,{T) = (2j[ + l)(2j 2 + 1) 

xk&^( T )e AE/kBT - (34) 

The ratio R Jlj2 ^(T) = k^^T)/^^ is 
proportional to an exponent with the argument depend- 
ing on AE. In this work we computed two direct-reverse 
processes in the HD + p-H 2 collisions, specifically: 

HD(0) + H 2 (2) ^HD(2) + H 2 (0), (35) 

with AE = 96.6 cm -1 . It would be useful to check how 
well the computed thermal rates ([8]) in this work satisfy 
the detailed balance equation (jMJ). Fig. 9 represents 
these results. It is shown that all results are in a satis- 
factory agreement with each other. 
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FIG. 9. (Color online) The ratio R Jlj2 ^j> j> (T) for the pro- 
cesses (|35|>. Computation with the use of both potentials: the 
BMKP and the modified DJ PESs. The open circles are the 
values of the exponential function from the right side of Eq. 



IV. CONCLUSION 

In this work a rotational method has been applied for 
the transformation of the 4-dimcnsional rigid monomer 
model H2-H2 PES of Ref. [H to the non-symmetrical po- 
tential appropriate for calculations of the HD+H2 col- 
lisions. Different low energy elastic and state-selected 
inelastic cross sections as well as the thermal rate coef- 
ficients for HD+H2 have been computed and compared 
with previous calculations, where available. The rota- 
tional energy transfer in HD+H2 is of importance for 
the thermodynamics of the ISM^. By now few and 
rather conflicting results are available for the low energy 
HD+H2 rotational energy transfer, see for example^ and 
references therein. The BMKP PES^i has been already 
applied to HD+H^^. However, this PES may have 
failures. The fact was mentioned in Refs. [TH and|l3| and 
ir*2i. Therefore in this paper a new attempt has been 
undertaken to carry out alternative computational meth- 
ods for HD+H 2 collision. In case of the BMKP and DJ 
PESs the necessary steps for each potential have been 
described in Sees. Ill Bl and III CI and also in Ref. H3- 
In the case of the BMKP potential, which is a full six- 
dimensional surface^, the transformation from H2-H2 to 
the H2-HD system was done by shifting the center of 
mass in one H2 molecule to the center of mass of the 
HD molecule. Because the DJ PES has been formulated 
for the rigid monomer rotor model, the transformation 
methodology was more complicated. Simply, the R\ and 
i?2 coordinates arc not available in this case, they have 
fixed lengths. In this paper the transformation has been 
accomplished by rotation of the space fixed OX YZ co- 
ordinate system, i.e. by the redirection of the R3 vec- 
tor. The new vector R' 3 connects the center of masses 
of the H2 and the HD molecules as shown in the Fig. 
f2J This procedure obtains new coordinate angles for the 



Jacobi vectors R±, R2, and R3, and a new PES as in Eq. 
(P0"|) . New experiments that measure the state-to-state 
rotational cross sections in the HD+0-/P-H2 collisions at 
low temperatures are needed. Thereafter theoreticians 
and astrophysicists would be able to compare compu- 
tational results with available experimental data. This 
type of work was recently accomplished, for instance, 
for the par C1-H2-I-H2 collision^. Here it would be use- 
ful to mention other contributions on hydrogen-hydrogen 
collisions^—. Additionally, with the use of new HD+H2 
results for the thermal rate coefficients one could carry 
out new computation of the HD-cooling function men- 
tioned in the introduction^. 

In conclusion, another interesting system worth men- 
tioning is HD+HD. For this collision there are relatively 
old experimental state-to-state rotational probabilities 
for a few selected states^. These old data can be useful 
in comparisons with the computational results obtained 
with different H4 PESs: such as the available DJ and 
the BMKP PESs or some relatively new potentials, for 
example from work o 14 i 52 . In the case of the DJ surface 
it would be possible to again apply the rotation proce- 
dure of the OXYZ coordinate system as performed in 
this paper. 
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